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Ab initio molecular dynamics simulation is used to study the structure and electronic properties 
of the liquid Ga-Se system at the three compositions Ga2Se, GaSe and Ga2Se3, and of the GaSe 
and Ga2Se3 crystals. The calculated equilibrium structure of GaSe crystal agrees well with available 
experimental data. The neutron-weighted liquid structure factors calculated from the simulations 
are in reasonable agreement with recent neutron diffraction measurements. Simulation results for 
the partial radial distribution functions show that the liquid structure is closely related to that of 
the crystals. A close similarity between solid and liquid is also found for the electronic density of 
states and charge density. The calculated electronic conductivity decreases strongly with increasing 
, Se content, in accord with experimental measurements. 
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I. INTRODUCTION 



The study of liquid semiconductor alloys goes back many years, and has been a splendid source of iusights into the 
relationships between atomic ordering, chemical bonding and electrical properties in condensed matteia. In some well 
studied systems like Cs-Auo and Mg-Bil3, where the elements have rather different electronegativities, the electrical 
' conductivity passes from highly metallic values to insulating values, as a gap at the Fermi level opens up in the 
" density of states and the bonding changes from metallic to ionic at the stoichiometric composition. In cases where 
one element is a metal and the other a semiconductor in the liquid state, the chemical bonding can be continuously 
tuned from metallic to covalent as the composition is varied. The Ga-Se system studied in this paper is an example 
of the latter kind. |-. 

We have used ab initio molecular dynamics simulatioru (AIMD) to study the liquid Ga-Se alloy at the three 
compositions Ga2Se, GaSe and Ga2Se3. The simulations have been used to study the liquid structure, the electronic 
density of states and the spatial distribution of electrons as a function of composition. We shall show that the 
relationship between the properties of the solid and the liquid is very instructive, and we have used the same methods 
to study the two stoichiometric solids GaSe and Ga2Se3. 

The liquid Ga-Se system has been rather little studied experimentally, and we chose to work on it mainly because 
we were aware of plans to perform neutron diffraction measurements on the systemQ. These measurements, as well as 
studies of the electrical conductivity, were completed while the work was in progress, and we shall present comparisons 
^ , with both sets of results. The system is also expected to showppuite a close resemblance to the liquid alloys In-Se (the 
• electronic properties have been, studied by Okada and OhnoO), Ga-Te (structure measurements have been done by 
Takeda et alB and Hoyer et alH^d electronic properties by Valiant and FaberB) and In-Te (structure measurements 
have been made by Hoyer et al. 113, conductivity and thermopower by Popp et al. O). 

It is useful to consider what structural and electronic changes would be expeejted as the composition of £-Ga-Se is 
varied. Liquid Ga is a typical liquid metal, with a coordination number of ^ 9tHI, while £-Se has a covalently bonded 
chain-like structure with a coordination number of 2 (Ref. |l^ . It is clear then that there must be radical changes of 
structure across the composition range. However, the electronegativity difference between Ga and Se is fairly small, 
so that strongly ionic bondingjs not expected. This is confirmed by the structure of the crystals. In the rather 
unusual layer structure of GaSellil, the Ga atoms arc tetijaJacdrally bonded to Se and to each other, and the Se atoms 
have threefold coordination to Ga. The Ga2Se3 crystalHl shows tetrahedral coordination of both elements. In the 
middle concentration range of the liquid, we might therefore expect a covalently bonded low-coordinated structure 
dominated by Ga-Ga and Ga-Se bonds. This expectation is fully confirmed by our simulations. 

The AIMD technique is ideally suited to this type of problem. The basic principle is that the total energy of the 
system and the forces on the atoms for any atomic arrangement are calculated by solving Schrodinger's equation for 
the valence electrons to find the self-consistent ground state. By simulating the system in thermal equilibrium, we 
therefore generate the liquid structure, the chemical bonding and the electronic structure in a completely unified way, 
and without any prior assumptions or adjustable parameters. In the past few years, AIMD has become a widely used 
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technique for studymg liquid metals and semiconductorstStZl. The work reported here is closely related to our recent 
simulations of ^-GaBand ^-Ag-Seta. 

The paper is organized as follows. In the following section, we summarize the AIMD techniques we use. Sec. 3 
presents our calculations on the equilibrium structure and electronic structure of the crystals. Our AIMD simulations 
on the liquid alloys are then reported in Sec. 4. The significance of the results and comparisons with related liquid 
alloys are discussed in Sec. 5. 



II. TECHNIQUES 

The first-principles simulation methods used here are the same as those used in our work on ^-GaEl, and we give 
only a brief summary of the main points. Only valence electrons are treated explicitly, and the core states are assumed 
to be identical to those in the free atoms. For the Ga-Se system, all states up to and including the 3d states in both 
species are treated as part of the core. The interactions between valence electrons and the cores are represented by 
norm-conserving nonlocal pseudopotentials, which are constructed ab initio via calculations on the free atoms (see 
below). The calculations are performed in periodic boundary conditions, with the electronic orbitals expanded in plane 
waves, all plane waves being included whose kinetic energy is less than a chosen cutoff i?cut- The exchange-correlation 
energji-is represented by the local density approximation (LDA), the form used here being that due to Ceperley and 
Aldeio. In the dynamical simulations, the ions follow classical trajectories determined by the Hellmann-Feynman 
forces obtained from the first-principles calculations, while the electronic sub-system remains in the ground state at 
every instant (the Born-Oppenheimer principle). 

The problems p£pciiforming accurate first-principles simulations on metallic systems have been extensively discussed 
in the literaturellaOiH3. We use here the Fermi-smearing technique with orbital occupation numbers treated as 
dynapiical variables. The smearing is handled with the quasi-Gaussian smearing function introduced in our ^-Ga 
workllS. Minimization of the (free) energy to obtain the self-consistent ground state for each ionic configuration is 
performed by the preconditioned conjugate gradients method. The calculations were performed with the all-bands 
version of the CASTEP codecJ running on the Fujitsu VPX 240 at Manchester. 

The norm-conserving pseudopotential fotnGa is identical to the one used in our ^-Ga work. For Se, we generated a 
pseudopotential using the standard Kerkercj method. The s and p components of the pseudopotential were generated 
using the atomic configuration 4s^4p'', and the d component using the configuration 4s^4p^-^^4c?"-^^. The core radii 
were taken to be 2.0, 2.0 and 2j3,a.u. for s, p and d components respectively. The pseudopotentials were represented 
in the Kleinman-Bylander formtj, witk.the p-wave treated as local for both Ga and Se, and the non-local parts of the 
pseudopotentials treated in real spaceEII. r— . 

The electronic densities of states (DOS) for the solids were calculated using the standard tetrahedron methodcS. 
For the liquids we used two approaches. The first method was to use F-point sampling only and to average over all 
ionic configurations, and the second was to use many fc-points for a few chosen configurations. Both methods gave 
essential the same DOS. 



III. CRYSTALLINE PHASES 

The /3-GaSe crystal has the layered structure shown in Fig. |^, with every layer consisting of four planes of atoms. 
Each of these planes contains one type of atom, the ordering of the planes being Se-Ga-Ga-Se. The symmetry of the 
crystal is hexagonal, with the hexagonal axis perpendicular to the layers, and the atoms in each plane are arranged in 
a regular close-packed hexagonal lattice. The interatomic bonding is strong within each layer, but weak between the 
layers, and this means that there are different ways of stackiag the layers which have almost the same energy. Four 
polytypes differing in the stacking are experimentally knownll^, and the observed form depends on the preparation 
method. There has been controversy about which polytype is most stable. 

Our aims in studying the crystal are to test our pseudopotentials and to obtain an understanding of the electronic 
structure, so that our main concern is with the strong intra-layer bonding. For this purpose, it does not matter which 
polytypCj-we examine, and we have chosen to work with the /9-phase, for which rather precise diffraction data are 
availablecj. Because of the relation between adjacent layers in this polytype, the primitive cell contains a total of eight 
atoms. The crystal structure is characterized by the two lattice parameters a and c and the two internal parameters 
u and V. The u and v parametecs specify the positions of the Ga and Se planes respectively along the ifeaxis. (For a 
detailed definition, see WyckoflEj; the z-coordinates for Ga and Se used in the paper of Benazeth et alSB correspond 
exactly to these u and v parameters.) 
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The calculations on /3-GaSe were performed using a plane-wave cutoff of 250 eV. This cutoff was chosen on the 
basis of systematic tests on the convergence of the total energy as a function of cutoff, which showed that thej^sidual 
error is roughly 1 meV per 8-atom cell. Brillouin-zone sampling was performed using 18 Chadi-Cohen point£j in the 
full zone. Again, tests done with different fc-point sets indicate that the residual error is no more than a few meV per 
cell. 

We have performed a full structural relaxation of the crystal, which gives the following equilibrium parameters 
(experimental values in parentheses): a = 3.64 A (3.750 A), c = 15.76 A (15.995 A), u = 0.174 (0.1736) and v = 0.602 
(0.6015). The reasonably close agreement with experiment fully confirms the adequacy of our pseudopotentials. A 
slight underestimation of lattice parameters is commonly found in DFT calculations, and we note that the calculated 
lattice parameter in our earlier work on crystalline Ga was also too low b^-p^ 3 %. During the relaxation process, we 
have noted the expected strong anisotropy in the energetics of the crystalEd, with changes of the c-parameter causing 
much smaller energy changes than those of the a-parameter. 

The electronic DOS calculated for the equilibrium structure is shown in Fig. |2| (upper panel) r— .The calculations 
indicate that the material is a semiconductor, with a band gap-pf 1.35 eV (experimentally about 2 eVE2l). The electronic 
structure of /3-GaSe was studied many years ago by Schliiteifij using an empirical model pseudopotential, and we find 
that his band structure is in semi-quantitative agreement with our results. So far as we know, ours are the first ab 
initio calculations on the electronic structure of GaSe. 

The valence part of the DOS consists of an isolated peak at - 14 eV (bands 1-4), a double peak at - 8 eV (bands 
5-6 and 7-8), and a broader distribution extending from - 6 eV up to the valence band maximum (VBM) (bands 
9-18). A clearer understanding of these features can be obtained by studying the partial electron densities associated 
with chosen energy ranges. These densities (Fig. ^ show that the peak at - 14 eV arises from Se(4s) states, while 
the double peak arises from bonding and anti-bonding states associated with Ga-Ga pairs, but with some weight on 
neighboring Se atoms. These states can be regarded as made of Ga(4s) states, with a strong admixture of other 
states. The broad peak in the DOS between - 6 eV and the VBM appears to be mainly responsible for bonding 
in the material, and consists of Ga(4p) and Se(4p) states. This general analysis is consistent with what was found 
by SchliiterEi rJt is interesting to note that there have been recent first-principles calculations on the closely related 
material InSeES, and the DOS calculated for that system shows all the same qualitative features as ours. 

We have also performed calculations on the Ga2Se3 crystal, though here the experimental situation is rather 
unsatisfactory. Again, there are a number of closely related structures, all of which appear to be based on a defective 
zinc-blende structure, in which one third of the cation sites are vacant. One of these (7-Ga2Se3, wiiiph is the high 
temperature phase) has full cubic symmetry, and the distribution of vacancies appears to be randomllj. 

Our Ga2Se3 calculations were performed on a periodic system having an 80-atom cell (Ga32Se48) consisting of 
twelve adjacent zinc-blende cubes, with sixteen, randomly chosen, vacant cation sites per cell. For this system we 
have calculated the electronic density of states (without relaxing the system) and the valence-charge density. The 
calculations have been performed with the same plane-wave cutoff as for /3-GaSe. 

The electronic DOS is shown in Fig. || (lower panel). Again, we are dealing with a semiconductor, the gap in this 
case being about 0.6 eV. The valence DOS consists of the Se(4s) peak at about -14 eV and a broader distribution 
from -10 eV up to the energy gap. The peaks in the DOS are not so well visible as in the case of /3-GaSe because of 
the randomness of 7-Ga2Se3. 

IV. LIQUID ALLOYS 

A. Structure and dynamics 

We have performed AIMD simulations of the liquid Ga-Se system for the three compositions Ga2Se, GaSe and 
Ga2Se3 at the temperature 1300 K. All the simulations were done on a repeating system of 60 atoms using F-point 
sampling, with a plane-wave cut-off of 150 cV. This cut-off is less than the one used for calculations on the solids, 
but is adequate for simulations of the liquids. The time step was taken to be 3 fs, and the Fermi smearing was 
0.2 eV. The density of these liquid mixtures is not experimentally known, and the densities we use are obtained from 
experimental values for £-Ga. and £-Se. Data for the density of i-Ga are available up to ~ 1000 K, and the value 
at 1300 K can be estimated by a small extrapolation. The case of i-Se is more tricky, since the boiling point under 
atmospheric pressure is only 958 K. We have estimated the (hypothetical) density at 1300 K by linear extrapolation 
of experimental data between 490 and 958 K. Finally, the densities of the liquid mixtures at 1300 K are obtained 
by linear interpolation between the estimated values for ^-Ga and ^-Se at this temperature. Some time after the 
work was staricd, a cross-check became possible against recent neutron diffraction data on ^-Ga-Se at the same three 
compositiongJ, and this showed that our estimated densities were correct to within better than 5 % in all cases. It 
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is worth mentioning that hnear interpolation between c-Ga and c-Se predicts the densities of sohd GaSe and Ga2Se3 
very accurately. 

We initiated the simulations by starting from a typical configuration of atoms taken from our previous £-Gsl 
simulationalj, and replacing some of the Ga atoms by Se atoms. At each composition, our simulated system was 
equilibrated for 1 ps, and then a production run of 4 ps was generated. The temperature and liquid structure were 
monitored throughout the simulations, and the stability of these quantities indicated that the 1 ps equilibration period 
was adequate. 

The structure of the simulated liquid can be compared directly with that of the real system through the static 
structure factors Sapik). These quantities, which give a measure of the intensity of density fluctuations as a function 
of wavevector fc, are defined by: 

^a/3(fc) = (P«(fe)p;(fe) ) . (1) 

Here the dynamical variable pa{k) representing the Fourier component of the atomic density of atoms type a at 
wavevector k is given by: 

p{k) = N-^/^ £ exp(zfe . r,) , (2) 

where is the position of atom i and Na is the number of atoms of type a in the system. The angular brackets 
in Eq. (|l|) denote the thermal average, which in practice is evaluated as the time average over the duration of the 
simulation. The structure factors were calculated for vectors k compatible with the periodic boundary conditions. In 
practical calculations of Sapik), we also average over k vectors having the same magnitude. The neutron weighted 
structure factor Sn{k) measured in a neutron diffraction experiment is given by: 

Sn{k) = -5 , (3) 

where Cq and ba are the concentration and the neutron scattering length of species a respectively. In the present 
work, we took the scattering lengths for Ga and Se to be 6Ga=7.288 fm and 6se=7.97 fm (Ref. |^). We shall also need 
to refer to the total structure factor S{k), which is defined by: 

S{k)^Y.V^S^t^ik) ■ (4) 

a/3 

In fact, S{k) is virtually identical to Sn{k) for the present systems, since &Ga and 6se are almost the same. 

The neutron- weighted liquicL^tructure factors Sn{k) calculated for the three compositions are compared with the 
recent neutron diffraction dataEI in Fig. ^. The agreement is good for wavevectors above ~ 3 A~^, and is reasonably 
satisfactory below that. Since the main discrepancies are at low fc, it is possible that the rather small size of our 
simulated system may be infiuencing the results; limitations of computer power prevent us from going to larger systems 
at present. The main changes in the experimental S'„(fc) with increasing Se content are the appearance of a small 
peak at ~ 1 A""*^, a shift and enhancement of the peak at ~ 2 A""'^ and an increase in amplitude and shift in phase 
of the oscillations beyond ~ 3 A~^. All these effects are correctly reproduced in the simulations. 

The origin of the features in 5'„(fc) can be understood by examining the partial structure factors Saf3{k) shown in 
Fig. ||. These results make it clear that the main peak in S'„(fc) at ~ 3 A"-"^ arises from a superposition of almost 
coincident peaks in all three Sap{k). The peak at ~ 2 A~^ comes mainly from a peak in Sse-So{k) , which is partially 
cancelled by a trough in S'Ga-Sc(fc); both of these features shift to lower k as the Se content is increased. The 
appearance of the small peak at ~ 1 A~^ is due mainly to a growth of Ga-Se correlations with increasing Se content. 

The real-space structure of the liquid can be understood through the partial radial distribution functions (RDFs) 
9ai3{r) (Fig. g), and the coordination numbers of their first peaks (table ^. One can notice very strong similarities 
between the structures of solids and the corresponding liquids. Not only are the interatomic distances very similar 
but also the coordination numbers are very close. 

It is useful to relate the forms of the RDFs to the corresponding crystal structures. The main peak in (?Ga-Se arises 
from directly bonded Ga-Se pairs, and its distance is close to the corresponding bond length in the crystals. The main 
peak in gso-Se comes from pairs of Se atoms bonded to the same Ga atom, and its distance remains unchanged for 
the same reason that the Se-Se distance is almost the same in the GaSe and Ga2Se3 crystals. We note, however, the 
formation of a small peak in gsc-Sc at 2.4 A at the Ga2Se3 composition, which is due to Se atoms dirjCctly bonded 
to each other; this is confirmed by the fact that 2.4 A is almost exactly the bond length in crystallineEHl and liquid 
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Sella. At low Se content, gca-Ga consists mainly of a peak at ^ 2.6 A, whose height decreases strongly on going to 
Ga2Se3, with the growth of a second peak at ~ 3.8 A. This is related to the presence of Ga-Ga bonds in crystalline 
GaSe and their absence in Ga2Se3. 

We have studied the diffusion behavior of Ga and Sc by calculating the time-dependent mean square displacements 
for the two species. These show the usual liquid-like behavior, with a rapid transition to a linear dependence on time 
after about 0.1 ps. The self-diffusion coefficients -Dca and D^c are obtained in the usual way from the slope, and the 
calculated values are shown in Table |l[ The results show that both diffusion coefficients decrease rather strongly with 
increasing Se content, presumably because of the formation of a faii% stable covalently bonded network. Interestingly, 
though, recent AIMD simulations that we have made on pure £-S£3 show that its diffusion coefficient at 1375 K has 
a high value of 1.5x10""* cm^s~^. We conjecture that Dgc has a minimum at the composition Ga2Se3 and increases 
rapidly thereafter. 

B. Electronic properties 

The total electronic densities of states (DOS) of the GaSe and Ga2Se3 liquids are compared with our results for 
the corresponding crystals in Fig. |^. The comparison shows the very close resemblance of the solid and liquid phases, 
which might be expected from the similarity of their short-range order. In both cases, the liquid-state DOS is simply 
a rather broadened version of the crystal DOS but shifted to higher energies. The main difference for GaSe is the 
replacement of the band gap by a minimum in the DOS. The variation of electronic structure with_jcomposition 
is illustrated in Fig. |8[ where we include the DOS for £-Ga at 982 K taken from our previous workll3. The ^-Ga 
system is highly metallic, and has a free-electron-like DOS. The change on going to Ga2Se is very marked, with the 
appearance of the feature due to Se(4s) states and the formation of the broad two-component distribution consisting 
of hybridized Ga(4s/p)-Se(4p) states characteristic also of GaSe and Ga2Se3. The DOS at the Fermi level shows a 
monotonic decrease with increasing Se content, reaching a zero value at the Ga2Se3 composition. 

The close resemblance between the electronic structure of the solid and liquid is confirmed by a study of the electron 
density distribution. A useful way to do this in the liquid is to examine the density on a plane passing through two 
neighboring Ga atoms and a Se atom neighboring one of these Ga atoms. We show in Fig. ^ the energy resolved 
partial densities for a typical configuration in £-GaSe, which can be directly compared with the corresponding results 
for the crystal (Fig. ||). All the characteristic features found in the solid - the Se(4s) states (bands 1-30), the bonding 
and anti-bonding states on the Ga pair (bands 31-45 and 46-60), and the Ga-Se and Ga-Ga bonding states associated 
with the upper part of the valence band (bands 61-135) - are clearly visible in the liquid. We find a similar close 
resemblance for solid and liquid Ga2Se3, as illustrated in Fig. |l^. 

The relation between the electron density distribution in the solid and liquid can be pursued to a more quantitative 
level. To do this, we have chosen particular Ga-Ga and Ga-Se pairs in liquids of different compositions and in the 
solids, and plotted the electronic density along the bonds. Since the distances between neighbors vary in the liquid, 
we reduce lengths to a common scale by dividing by the interatomic distance. Some typical comparisons are shown in 
Fig. |ll|. The agreement between the solid-state and liquid-state curves, even for different compositions, is remarkable, 
and indicates that there is almost no change in the short-range electronic structure between the different phases. 

We haj^e calculated the electronic d.c. conductivity of the three liquid phases ijUg the Kubo-Greenwood approx- 
imationO. The technique used is the one described in several previous paperalj'tj. The calculations were done by 
averaging over the full duration of each simulation, and used F-point sampling, which is expected to be adequate 
for present purposes. Fermi-Dirac thermal occupation numbers were included in. the calculations. The results are 
compared with the very recmit experimental measurements of Lague and BarnesQ in Fig. |l^. In this graph we also 
include results for pure Galla. Given the approximations involved, the agreement is as good as can be expected. 
The dramatic decrease of conductivity is, of course, associated with the reduction of DOS at the Fermi level already 
described. 



V. DISCUSSION AND CONCLUSIONS 

The good agreement between the calculated equilibrium structure of the GaSe crystal with available data, as well 
as the satisfactory comparisons of the liquid structure factors with recent diffraction data, give confidence in the 
realism of our AIMD simulations. One of the main findings from these simulations is the very close relation between 
the properties of the solids and the liquids. The positions of the peaks in the liquid RDFs for GaSe and Ga2Se3 
are close to the interatomic distances in the crystals, and the coordination numbers are also similar. As suggested 
in the Introduction, the bonding in the liquids appears to be a mixture of metallic and covalent, with the structure 
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dominated by Ga-Ga and Ga-Se bonding in the composition range we have examined. These findings about the 
similarity of sohd and hquid and the. nature of the bonding are in hne with indications from early diffraction work on 
the related Ga-Te and In-Te liquidaJ. 

Clearly at some point in the composition range, direct bonding between Se atoms must become important, since 
the structure of pure ^-Se is entirely determined by this bonding. According to our simulations, this direct bonding 
is insignificant in the range up Ga2Se3. However, there is clear evidence in the gsc-Se RDF for the beginning of 
this effect at the. Ga2Se3 composition. It is interesting here to compare with our recent AIMD simulations on the 
£-Ag-Se systenO, where we found that direct Se-Se bonding begins at the Ag2Se composition, and rapidly becomes 
a dominant effect for higher Se contents. This strongly suggests that the point at which Se-Se bonding begins is 
determined by the maximum valency of the other component. In this context, ab initio simulations or diffraction 
measurements on systems such as Cu-Se and Cd-Se would be extremely interesting. 

The strong similarity of solid and liquid is also very clear from our results for the electronic DOS, where we have 
shown that for both GaSe and Ga2Se3 the Hquid DOS is essentially a broadened version of the soUd DOS, except for 
the disappearance of the band gap in the case of GaSe. Interestingly, the band gap does not disappear for ^-Ga2Se3. 
This behavior of the DOS is closely related to the composition dependence of the electrical conductivity, which, 
according to our rather limited results, decreases rapidly and monotonically as one passes from £-Ga to £-Ga2Se3, as 
is also found experimentally. 
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FIG. 1. Structure of the /3-GaSe crystal. Dark and light spheres represent Ga and Se respectively. 

FIG. 2. Calculated electronic DOS for /'3-GaSc (upper panel) and 7-Ga2Se3 (lower panel). For presentation purposes, the 
calculated DOS are convoluted with a Gaussian of width 0.1 eV. The horizontal scale represeritsthe difference of the energy E 
and the Fermi energy Ef. 

FIG. 3. The band-resolved and total charge densities for the /3-GaSe crystal (units; 10^^ A^'*). Ga and So sites are marked 
by filled circles and squares respectively. Lengths are marked in A units. The first four panels show bands resolved charge 
densities corresponding to features in the DOS (for meaning of bands see text). The last panel shows the total charge density. 

FIG. 4. Simulated (solid lino) and experimental (circles) neutron-weighted structure factors of liquid Ga-Se at three compo- 
sitions. The vertical scale refers to the GaSe results, with Ga2Se and Ga2Se3 results shifted respectively up and down by one 
unit. 

FIG. 5. Faber-Ziman partial structure factors Sag{k) and total structure factor (see text) obtained from simulations of liquid 
Ga-Se at three compositions. 

FIG. 6. Partial radial distribution functions Qagir) obtained from simulations of liquid Ga-Se at three concentrations. Solid 
and dashed vertical lines indicate interatomic distances in the /3-GaSe and 7-Ga2Se3 crystals respectively. 

FIG. 7. Comparison of the electronic DOS for solid ( ) and liquid ( ) GaSe and Ga2Se3. 

FIG. 8. Comparison of the electronic DOS calculated for l-Ga, and £-Ga-Se at three compositions. 

FIG. 9. Band-resolved and total electron densities for i?-GaSe (units: 10^^ A^"^). Ga and So sites are marked by filled 
circles and squares respectively. Lengths are marked in A units. The first four panels show bands resolved charge densities 
corresponding to features in the DOS (for meaning of bands see text). The last panel shows the total charge density. 

FIG. 10. Total electron density for solid (left panel) and liquid (right panel) Ga2Se3 (units: 10"^ A~^). Lengths are marked 
in A units. 

FIG. 11. Electron density (units: 10"'^ A~"^) along lines joining neighboring Ga-Ga and Ga-Se pairs in various solid and 
liquid phases. 
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FIG. 12. Comparison of calculated (open circles) and experimental (filled diamonds) electronic conductivities for liquid 
Gai-xScx alloys. 



TABLE I. The average coordination numbers for liquid Ga-Se alloys and positions of the first peaks in the radial distribution 
fuction, in brackets the interatomic distances and the coordination numbers of the corresponding solid phase. Distances in A. 



liquid rCa-Ga WCa-Ga roa-Se WQa-Se 


WSe-Ga rse-Se 


nSe-Se 


Ga2Se 2.54 2.41 2.44 1.46 

GaSe 2.51 (2.44) 1.12 (1) 2.44 (2.45) 2.66 (3) 
Ga2So3 2.46 (3.85") 0.62 (0) 2.43 (2.35) 3.39 (4) 


2.92 3.95" 
2.66 (3) 3.95" (4.75") 
2.26 (2.67) 3.95" (3.85") 


6.2" 
8.0" (9") 
9.0 (12 ) 


" these numbers correspond to the second coordination shell. 






TABLE II. Calculated diffusion coefficients (in 10" 


"^cm^s"^) for Ga-Se liquids. 




liquid Doa 


Dse 




GaaSe 9.8 
GaSe 4.6 
Ga2Se3 2.2 
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